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Introduction 


This investigation is concerned with the dynamic characteristics of the 
Long-Term Zonal Earth Energy Budget Experiment (LZEEBE) spacecraft (see Fiq. 1) 
The spacecraft mass center moves in a known circular orbit about the earth. On 
the other hand, the spacecraft attitude motion is expected to be as random as 
possible. Hence, capture in a gravity-gradient configuration must be avoided. 

The mathematical model consists of a rigid hub and three flexible booms 
equally spaced in a plane, referred to as the spacecraft plane, where the 
angle between any pair of booms is 120°. At the tips of the booms there are 
spherical balloons coated with different materials, so that the solar ra- 
diation produces not only forces but also torques. 

Assuming that the orbital motion is known, the dynamical problem reduces 
to that for the rotational motion of the spacecraft as a whole and the elastic 
motion of the flexible booms. Representing the displacements of the booms 
by the first bending mode of the associated fixed-base booms with balloons 
at their tips, the behavior of the spacecraft can be simulated by means of 
a nine-degree-of-freedom nonlinear system subjected to external excitation. 

No closed-form solution can be obtained for such a system, so that the equations 
of motion must be integrated numerically. 

The nine-degree-of-freedom simulation has been cast in a form suitable 

4 

for numerical integration and programmed for digital compuation. Computer 
results indicate that the spacecraft will perform as desired, provided the 
injection in orbit does not create conditions favorable to capture. 
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The Generalized Coordinates 


Let us consider an inertial system XYZ with the origin at the sun's 
center S; axes X and Y are in the ecliptic plane, with axis X along the 
vernal equinox, and axis Z is normal to the ecliptic plane. The earth's 
center E moves in the ecliptic plane in a circular orbit around the sun, 
so that at any time the position of E is defined by the radius vector 
making an angle A ^ with respect to the vernal equinox (see Fig. 2). The 
equatorial plane of the earth intersects the ecliptic plane along an axis 
X' parallel to X. Denoting by Y' the axis normal to X' and in the equatorial 
plane and by V the polar axis of the earth, and letting a be the angle 
of inclination of the polar axis relative to the ecliptic, we conclude that 
the relation between the systems X'Y'Z' and XYZ can be written in the 
matrix form 



where 

1 0 0 
0 COSa -sina 
0 sina COSa 



(1) 


( 2 ) 


plays the role of a rotation matrix (see Ref. 1). 

Next, we wish to define the orbit of the spacecraft. To this end, we 

c 

introduce the set of axes X"Y"Z" obtained from X'Y'Z' by a rotation A^ 
about axis Z'; axis X" defines the ascending node of the spacecraft orbit 
(see Fig. 3 ). The relation between X"Y"Z" and X'Y'Z' is simply 




X 

Y 

Z 




( 3 ) 
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where 



cos A^ sinA^ 0 
-sinA^ cosAj^ 0 
0 0 1 


(4) 


Denoting by i the orbit inclination relative to the equatorial plane, the 
spacecraft orbit plane can be defined by axes obtained from axes 
X"Y"Z" by means of the rotation i about X". The relation between the two 
sets of axes is 


H 

r n 

fx") 

Y"' 

= R(D 

Y" 

Z" 1 

L J 

Z" 

v ) 


/ 


where 



1 0 0 

0 cos i sin i 

0 -sin i cos i 


(5) 


( 6 ) 


The spacecraft moves in a circular orbit about E in the X"'Y'" plane. 
The position in the orbit of the spacecraft mass center C is defined by 
the radius vector from E to C, where the direction of the vector is 
defined by the angle <p = fit measured from the ascending node. It will 
prove convenient to introduce a set of "orbital axes" abc with the origin 
at C and with axis a along Rq, axis b tangent to the orbit and in the di- 
rection of motion, and axis c normal to the orbit plane. The relation be- 
tween systems abc and x'"Y'"Z " 1 is simply 




where 


r , vi 

COSij; 

si nip 

0 

RU) = 

-sin^ 

COSiJi 

0 

L J 

0 

0 

1 


(7) 

(8) 
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The overall motion of the spacecraft can be conveniently described by 
a set of body axes xyz, defining the rotational motion of the spacecraft 
relative to axes abc, and by the elastic deformation of any point on the 
spacecraft relative to xyz. Assuming that the triad xyz, initially coincident 
with axes abc, undergoes the rotations e 2 about y, -e^ about x, and e 3 
about z in that order (see Fig. A), the relation between xyz and abc can 
be shown to be 

( x) (a 

y •=[*]< b . (9) 

' ' i 

where [«,] is the matrix of the direction cosines between the two sets of 
axes; its explicit form is 

ce 2 ce 3 - se.| se 2 se 3 ce-j se 3 - se 2 ce 3 - se-j ce 2 se 3 

- ce 2 se 3 - s e i se 2 ce 3 ce^ ce 3 se 2 se 3 - se-| ce 2 ce 3 

ce-j se 2 se i ce-| ce 2 

where se • = sine, and ce. = cose. (j=l,2,3). Note that of all the quantities, 

J J J J 

introduced to this point, only e^(t), e 2 (t), and e 3 (t) are generalized 
coordinates, as all the remaining quantities are either constant or they 
are known functions of time. 

It follows from the above that the relation between the body axes xyz 
and the inertial axes XYZ can be written in the compact form 

(ID 
( 12 ) 


X 

= CL] Y 
Z 

where 

[L] = [i][R(*)][R(1)][RU N )][R(a)] 
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is an overall matrix of direction cosines. 

The elastic displacements are assumed to be caused by flexure and to 
take place in two orthogonal directions. Letting x^(i=l,2,3) be the nominal 
longitudinal directions of the booms we shall denote the elastic displacements 
in the plane of the spacecraft by v.(x. ,t) and those out of the plane by 
w i (x i ,t) (see Fig. 5). 


Lagrange's Equations of Motion in General Form 

The motion of the spacecraft can be adequately described by nine generalized 
coordinates, three rotations e.(t) (j=l,2,3) and six elastic displacements 

V/ 

Vjfx.'t) and w^x^.t) (i=l,2,3), where the elastic displacements depend not 
only on time but also on spatial position. For simplicity, it is assumed that 
the spacecraft mass center coincides with its geometric center at all times, i.e^ 
its position does not shift relative to the spacecraft because of elastic 
displacements. 

To derive Lagrange's equations, it is necessary to produce expressions 
for the kinetic energy, potential energy, and work function. First, we must 
define certain vector quantites. Denoting by r Q the position vector of any 
point in the hub relative to C, the position of that point relative to the 
inertial space is simply 


Bo 


Be + B c + c 0 


(13a) 


On the other hand, denoting the nominal position of a point on boom i relative 
to C by r. and the elastic displacement of that point by u^ , the position of 
the point relative to the inertial space is 


5i = h + Sc + ri + 


i = 1,2,3 


(13b) 
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For convenience, let us assume that and are expressed in terms of com- 
ponents along the inertial system XYZ and that r Q , r^ , and u^ are in terms 
of components along the body axes xyz. If the body axes xyz rotate with 
the angular velocity o> relative to an inertial space, the absolute velocities 
of the point in question are 


So 


S E + Sc + 


“ * r 0 


(14a) 


and 


= Be + ~C + yj + 'to x (r.j + u^) , i - 1 , 2,3 


(14b) 


where u- denotes the elastic velocity of any point on boom i relative to 
axes xyz. 

The kinetic energy can be written in the general form 
3 


1 = 1 I 

L 1=0 'ml 


R . *R-f dm^ 


(15) 


Inserting Eqs. (14) into (15), we obtain 

3 r 

T = i- m(R r + R r ) * (R r + R r ) + \ £ 

2 -E ~C ~E ~C 2 i=Q J m , 

1 


[to x (r i + u. )]*[aj x (r i + u i )]dm i 


(16) 


m,- 


(Ci 


+ (r. + u.) x ul-dm. + E u_!*u! dm. 

~ i=l 'mv 1 i=l 'm\ 1 1 1 

3 3 

where m = E i=0 m^ is the total mass of the spacecraft. Moreover, e.._q 

+ ._ Uv)^.dm. - ’= -0C by the assumption that the center of mass of the spacecraft 

does not shift relative to the body axes xyz during motion. Introducing 

the notation 


3 

K = E 
: i=i i 


"K 


(Ti + y-j) x yj 


(17) 


where K is recognized as the angular momentum about C due to elastic velocities 
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alone, and letting J be the inertia dyadic of the spacecraft in deformed 
state about xyz, the kinetic energy can be written in the compact form 


T ■ T C + T rel 


where 


T c 2 ^(?£ Bq) * ^~E ^ B C } 


( 18 ) 


(19) 


is the kinetic energy due to the motion of C and 

3 


T rel = + + 


m. 

i 


u:-u! dm.j 


(20) 


is the kinetic energy due to the motion of the spacecraft relative to C. 

Because Tq contains no generalized coordinates or velocities, it will be 
dropped in future discussions. Consequently, ‘the subscript rel will be 
dropped in Eq. (20). Equation (20) can also be written in matrix form. In- 
deed, denoting the inertia matrix representing J by [J] and the column matrices 
associated with the vectors ay, K, and u.j by lay’}, {K}, and {u.1}, respectively, 
Eq. (20) becomes 


3 

T = i {uy} T [J]{oy} c.+.a {K} T {oy} +,Z f (u 4 > T {u: > dm, (21) 

* i=l ^m i 1 1 1 

At this point it appears desirable to specify some of the quantities 
in T. It should be pointed out that all these quantities are in terms of 
components about the body axes xyz. Letting i, j, and k be unit vectors 
along xy y, and z, respectively, and recalling that there is a 120° angle 
between any pair of booms, we can write the position vectors 


Ti 


x l~ 


~2 = “ 
C3 


/3 


T x 2~ + T x 2~ 


2 x 3l 


/3 . 

t x 3 j 


( 22 ) 
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Moreover, the elastic displacements are as follows: 


u-, = v-, j + w-|k 


u» = 


-3 


- 4 V - \ v 2 j + w 2~ 

4 v 3 - 'I v 3 i + w 3 !S 


( 23 ) 


so that the relative elastic velocities are simply 


u] = + ^k 

“2 = - 4 V “ \ V- + v 2 k - 

“3 = 4 v 3! - \ v 3i + w 3~ 


(24) 


Considering Figs. 2-5, the angular velocity vector can be written in the form 


ij) = <i)' x i + oiyj + to 2 k 


(25) 


where, assuming that = 0, the components of w are 


to x = -q(s0 2 c© 3 + se-j c© 2 S0 3 ) - 0-j c©^ + 0 2 cc0| c0 3 

w = fx ( s© 2 s0 3 - se-j c0 2 c0 3 ) + 0-j se 3 + 0 2 c0-| s© 3 (26) 

w z = nc0-j c0 2 + 0 2 s e 1 + e 3 

The moments and products of inertia of the deformed spacecraft are 


J = J + 
xx xxo 


nv 


(v-, 2 + W^ 2 )dm 1 + ^ [^/3 x 2 - v 2 ) 2 + w 2 2 ]dm 2 

■'mJo 


yy 


+ x 3 + v 3^ 2 + w 3 2 ^ dm 3 

Jm -3 

l m , ( X 1 2 + w l 2 > dm , + J [*<*2 + + 

J t^( x 3 “ S 3 v 3 ) 2 + w 3 2 ]dm 3 


J + 

yy 0 
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J = J + 
zz zzo 


I (x-, 2 + v 1 2 )dm 1 + ^ f [(x 2 + ft v/ + {ft x 2 - v 2 ) 2 ]dm 2 


"1 "’2 

+ J-' C(x 3 - ft v 3 ) 2 + {ft x 3 + v 3 ) 2 ]dm 3 
% 3 

J xy = J xyo + | m : X T V 1^1 ' T m ^ x 2 + x 2 “ v 2 }dm 2 


+ 1 L {x 3 “ ^ V 3 )(/J x 3 + v 3 )dm 3 
lm 3 


J = J + 
xz xzo 


, x l w l dm l -?L (X 2 + /Jv 2 )w : 

mi ^ m« 


2 dm 2 " ? ^ 


m. 


(27) 

(x 3 - ft v 3 )w 3 dm 3 


J « J + , 

yz yzo j 


Vi dm i + ? 

m l 


(v? x 2 - v 2 )w 2 dm 2 - J J (ftx 3 + v 3 )w 3 dm 3 

m? 2 N 


where J xx0 » J yy 0 » ^y 20 are moments and products of inertia of the 
hub. Moreover, the elements of {K} have the expressions 


K x = 


m. 


« * 1 * * 

(v^w^ - w^v^dm-i + ^ [(^ x 2 - v 2 )w 2 - w 2 v 2 ]dm 2 


J m ?2 


- £ [(^ x 3 + v 3 )w 3 - w 3 v 3 ]dm 3 


K y = 


mi 


m 


x-jW^dm^ + £ 




• » 

[(X 2 + ft v 2 )w 2 - ft w 2 v 2 ]dm 2 


\ [ „ [(x 3 - ft V 3 )w 3 + ft WjCjJdBlj 
; mA 


(28) 


K z = x 1 v 1 dm 1 + 

; m-j 


x v 2 d ™ 2 + 

m 2 


, x 3 v 3 dm 3 
m 3 


The potential energy can be written in the general form 


V * V EL + V G 


(29) 
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where is the elastic potential energy having the expression (see Ref. 3) 

+ Ui t'( Ei -(^f) 2 f -^©) 2 ] 

+ ^ O 2 [' - 1(5) 2 ]h m 


Note that the first integral in Eq. (30) is due to axial forces and the 
second is due to flexure. On the other hand, Vg is the gravitational poten- 
tial energy, which has the matrix form (see Ref. 3) 

fi 2 T 

Vg - - — jf (t£ [J] - 3{& a } [J]{£ a )) (31) 

where (jO is the column matrix of direction cosines between the direction 

a 

of the vector Rg (which coincides with the direction of axis a) and axes xyz. 
Note that a constant term has been ignored in Vg. The matrix ’{«,' } has the 
explicit form 


<v ■ 


£ xa 

| 

C\J 

CD 

O 

i 

CO 

CD 

u 

CM 

CD 

to 

CD 

to 

00 

CD 

to 

l y a 

II 

o 

CD 

no 

1 

CO 

CD 

00 

to 

CD 

to 

CD 

PO 

CO 

CD 

o 

*za 

1 


r— 

CD 

O 

CM 

CD 

to 

J 


which is merely the first column of Eq. (10). 

It will prove convenient to eliminate the spatial dependence from v^ and 
w.j (i=l,2,3), i.e., to discretize the system. To this end, we use the assumed- 
modes method (see Ref. 2) and introduce the notation 

v-|(x-|,t) = 4>i (*i )q-| (t), w-j (x-j ,t) = <J»i (x*| ( t) 

v 2 (x 2 ,t) = <f> 2 (x 2 )q 3 (t), w 2 (x 2 ,t) = <j> 2 (x 2 )q 4 (t) (33) 

v 3 (x 3 ,t) = <{> 3 (x 3 )q 5 (t), w 3 (x 3 ,t) = 4> 3 (x 3 )q g ( t) 
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where $. (x i ) (i=l,2,3) are "shape functions". They can be taken as the first 
bending modes of the booms (see Ref. 2). The functions can be normalized 
so as to satisfy 


'm 


cj).. dm.j = 1 


1=1 ,2,3 


(34) 


Moreover, we can introduce the notation 


nr. 


x 1 * 1 dm i = b i , 


i=l »2,3 


(35) 


and 


f x? dm i = J i , i=l ,2,3 

•'mi 


(36) 


where is recognized as the moment of inertia of boom i about the z axis 
In view of Eqs. (33)-(36), Eqs. (27), (28), and (30) reduce to 

J xx = J xxo + l (J 2 i:J 3 ) ‘ 4 (b 2 q 3‘ b 3 q 5 ) + q l + q 2 + ^3 + q 4 + q 5 + ^6 
J yy = J yyo + ^ 4J i +J 2 +> V + ^ b 2 q 3' b 3 q 5^ + q 2 + 4*^3 + q 4 + 3^5 + q 6 

J zz = J zzo + J i + J 2 + J 3 + q l + q 3 + q 5 


J xy “ J xyo - + b l q l " ^2% ' l b 3 q 5 + “T^ q 3 “ q 5^ 


(37) 


J xz " J xzo + b l q 2 ‘ ^ b 2 q 4 + b 3 q 6^ “ 2^ q 3 q 4 “ q 5 q 6^ 
J yz = J yzo + % b 2% ~ b 3 q 6 } + q l q 2 ' ?< q 3 q 4 + q 5 q 6> 


• • i • • * • 

= ( c l2 ” ^2^1 ^ + 2 ^ 4^3 + ^2 ” ^3)^4 ” ^ 6^5 " ^3 ~ ^5)^5] 

K y = “ b l q 2 + 2^*^ q 4 q 3 + ( b 2 + ^ q 3^ q 4 + ^ q 6 q 5 + ( b 3 ' ^ q 5^ q 6^ (38) 

K z = b l q l + b 2 q 3 + b 3 q 5 
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and 



where it. was assumed- that EI^= , EI^ El ^ . (i=l, 2,3) i The. axial forces are 
due <to3 ceiitri.fuga.lt and ;gr&vitat16nalv£effects and can be obtained from Ref. 3. 
'Lagrange’s equations of motion can be written in the general form 



where e^. and are nonconservative generalized forces. In this particular 
case, they arise from solar radiation and internal damping. Equations (40) 
are second-order nonlinear differential equations and they are not very con- 
venient for integration purposes; it is actually more convenient to work 
with first-order differential equations. Our efforts will be directed next 
toward producing such a set of equations. 
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First-Order Differential Equations 


Considering Eq. (21) and recognizing that £ [ {ul} T {ul}dm. = 

6 2 i=l Jmi 1 1 1 

£ q*. , the second set of Eqs. (40) can be written in the form 

i=l 1 


3F (^)J' W> + {^7} T{ " } + q i ' { “ } 


9K ) T r i 3 V _ n 0 r 

‘ 3q.) 3q i ^i * 1 


(41) 


To transform Eqs. (41) into a set of twelve first-order equations, let us 
introduce the auxiliary variables 

P.j - q.j » i — 1,2,. ..,6 (42) 

In addition, let us introduce the 3x6 matrix [K*] defined as 


[K*] = 


9 ^x 9K x 


3K;. 


3 4; 9 q 0 


•1 




& ... 
aq-i 3q 2 3q 6 

3K^ 3K z . 3K Z 


[aq-i 3q 2 aq 6 

as well as the modified generalized forces 

3 K )T 


«r * q i - (at (%) 


3q. 


{«} 


+ 1<«} T 

2 3q i 


m 


_3 V 

9C H 


(43) 


(44) 


i=l;2,... ,6 

so that, arranging and Q* in column matrices, we can write Eqs. (41) and 
(42) in the form of the first-order differential equations 


{p} + [K*] T {«> =' {Q*} _ 


(45a) 


(q> = (pi 
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Next, let us turn our attention to the first set of Eqs. (40). The 
equations are in terms of the actual angular coordinates e,. It will prove 
more convenient, however, to work with a set of equations in terms of quasi - 
coordinates (Ref. 1). These equations can be written in the form 

(ar(lr)) + w'fS) *' + {N} r < 46 > 

where [<o/] is the skew symmetric matrix of angular velocity components, i.e. 



o 

-(O 

Z 

to 

y 

= 

<0 

Z 

0 

-to 

X 


-O) 

;y 

to 

"X 

» 

o 


(47) 


and {N}q and {N} r are the gravitational and radiation pressure torques, res- 
pectively. The gravitational torque can be shown to have the components 
(see Ref. 1 , p. 437) 


N x = 3Jr[(J -0 H sl + J s, i 
xg LV zz yy' ya za xy xa za 


N yG~ ^ J xx" J zz^xa A za + J yz £ xa £ ya 


N zG - 3ft [( J yy - J xx H xa * ya + ^xa^ya^za 


^xz^xaVa + J yz^za -2, ya^ 
J xy £, ya £ za + J xz^xa “ ^za^ 
J yz*xa £ za + J xy^ya “ £ xa^ 


(48) 


where a. , a, wa , are given by Eq. (32). The radiation pressure torque 
Xa ya Za 

is derived in the next section. 

From Eq. (21), we conclude that 


IS}- MM + fK> 

so that, introducing Eq. (49) into Eq. (46), we obtain 

+ [ J ] {<o } + { KJ + [o».][J]{(o} + Co,' ] {K> = {'N> g + (N} r 


(49) 


(50) 
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Next, we observe that 


{K} = [K*]{p> 


so that, introducing the notation 


{N*} = {N} + {N} r - [j]{o)} - [<o]([J]{u>} + {K}) (5J 

G K 

Eq. (50) reduces to 

[K*]{p}/+ = {N*} (5; 

On the other hand, from Eqs. (26), we conclude that the angular velocity 
vector can be written in the matrix form 


where 


{w} = [e]{e> + hu } 


-c6g ce^ s© 3 0 
[e] = se 3 ce^ c 0 3 0 
o se, 1 


Hse 2 ce 3 + se-j c© 2 se 3 )J 

U } s < se 2 se 3 - s e i ce 2 ce 3 > (56) 

/ C0 ^ c© 2 ■ y 

where {i } is recognized as the matrix of direction cosines between axis c 
and the body axes xyrz (see the last column of Eq. (10)). Equation (54) can 
be rewritten in the form 

{0} = [0] _1 ({io} - Q{£ c )) (57) 

Equations (53) and (57) represents another two sets of first-order differential 
equations replacing the set (46). 

Next, let us introduce the state vector 
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(58) 


{x} =■ 


i£ 

)l9> 

Tel 


e- 


which has the dimension eighteen. Then introducing the square matrix 


[A] = 


[l]‘; 

[K*f 

1 — 1 

0 

1 1 

i i — i 

1 o 
1 1 — 1 1 

[ k*]’T 

J. 

[J] 

! CO] ; 

i r — i i 
l O i 
i i i 

[0] : 

[0] 


i i — i i 
1 O i 
i i i i 

_[o] : 

[0] 

! to] : 

i — i 1 
j i 


as well as the column matrix of "generalized forces" 

( (Q*> 

{N*} 


{X} = 


Eqs. (45), (53), and (57) can be written in the compact form 
[A]{x> = {X} 


(59) 


(60) 


(61) 


where the elements of {X} are generally nonlinear functions of the state 
vector. Equation (61) can be rewritten as 


{x> = [A]”Vx> 


(62) 


which represents a form suitable for numerical integration. Note that 
Eq. (62) implies the existence of [A]" 1 , which in turn implies the existence 
of [e]~^ , with the further implication that [©] is nonsingular. From Eq. 
(55), we conclude that d.et [e] 0 provided cos 1 0. Hence, at values 

of e-j that are odd integer multiples of tt/2 Eq. (62) cannot be used, so that 
such values must be avoided in a numerical solution. 
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Radiation Forces and Torques 


There are several sources of radiation that can cause forces and possible 
torques on the spacecraft. These sources are: 

1. Direct solar photon radiation 

2. Solar radiation reflected by the earth and its atmosphere 

3. Radiation from the earth and its atmosphere 

4. Radiation from the spacecraft. 

The last three are one order of magnitude smaller than the first, so that they 
will be ignored. Moreover, it will be assumed that the solar photon radiation 
on the booms is negligible. Hence, the only forces and torques to be con- 
sidered are caused by solar photon radiation on the balloons. 

Next, we wish to obtain an expression for the force vector on a typical 
balloon. To this end, let us consider a differential element of area dA (see 
Fig. §) and denote by n the unit vector normal to the surface and pointing 
outward and by s the unit vector tangent to the surface and in the direction 
shown. Then, letting 0 be the angle between the incident radiation and n, 
the photon radiation force corresponding to the area dA is (see Ref. 4) 

dF = [-[(l+c rs )cose + §c rc j]n + (l-c rs )sine s]cose dA (63) 

o 

where I = energy per unit time through a unit cross sectional -area (in W/m ) 
c = speed of light 

c rs = coe ^ 1cier, t specular reflection 
c rc j = coefficient of diffuse reflection 
To obtain the force exerted on the balloon,; we must integrate Eq. (63) over 
the area, which is assumed to be spherical. From Fig. 6, we can write 
n = sine cos<j> 1^ + sine sin<j> 1 • + cose 1 

(64) 

s = cose cos<f 1_ + cose sin* 1 - sine V 

~K ~n 
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whpy»p 111 

~V -n 9 are unit vectors along axes £, n, C , respectively. More- 
over, the differential element of area has the expression 

dA = r^sine de d<f> (65) 

where r is the radius of the balloon. Denoting the force vector by 


F = FI + F 1 + FI 

n~n 


( 66 ) 


and integrating over the surface of the balloon, the force components can 
be shown to be 


F 


5 


F 


C 


= F =0 
n 

- o Jj- 1 i 2 \ 

- - 2 irr + 9 c rd ) 


(67) 


For further reference, we wish to calculate the forces on the various 
balloons in terms of components along the spacecraft body axes xyz. In view 
of Eqs. (67), we conclude that the force on any balloon has the magnitude 



_ p .,.21 /I 

- F ic " 2irr i c ( 2 


+ 9 c rdi^ 


1*1 .2.3 


( 68 ) 


where allowance has been made for the possibility that the balloons are of 
different sizes and that the coefficients of diffuse reflection are different. 
The direction of the force vectors coincides with the direction of the solar 
radiation. Hence, the direction of F. is parallel to the vector + R^. 
Because |R^| >> |R^| , it will be assumed that the direction is parallel to 
R^. This direction can be expressed in terms of inertial components as follows 



cos I + sin A^ J 


(69) 


where I and J are unit vectors along the inertial axes X and Y, respectively 
(see Fig. 2). It follows that 
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(70) 


El ■ 2 ’ r ^ C ( ? + f- c rdi )(cosi E ! + s,nA E J) 


The same vector can be expressed in terms of components along the spacecraft 
body axes xyz by writing 


El - F i* 1 + F iy i + F iz <£ 


(71) 


where the components F- x , F. , and can be obtained from Eq. (70) by means 
of the coordinate transformation (11). 

Next, we wish to use Eq. (71) and determine the radiation force and 
torque vectors {Q}^ and {Nl^. Ignoring the radius of the balloons compared 
to the length of the booms and using Eqs. (22), (23), and (33), we can write 
the position of the balloons as follows: 

r i ( -^ i ) (^ ft) - &i i 4>i Ui )<1-| (t) j ( ^■j )q2 ( t) k 

~ 2 ^ 2 ^ ^ H 2 () l 2 ’ t) - ~ 2 " ^2 ^ v/ ^" ^ 2 ^ £ 2 ^ 3 ^^~ 2 " ^2 ” ^ 2 ^ 2 ^ 3 ^^^ 

+ <t>2^ Jl 2^ c, 4^ t ^~ ( 72 ) 

r 3 (%) u 3 (& 3> t) — - 2" C^ 3 ~ ^3^3^5^^1 “ 2 " ^3 


+ <j> 3 U 3 )q 6 (t)k 


But the virtual work associated with the radiation forces can be written 
in terms of both actual coordinates and forces and generalized coordinates 
and forces in the form 


_ 3 

6W r = E F. • 6U.(£.,t) 
R i=1 ~1 1 1 


6 

i=l «1R «1 


(73) 


Inserting Eqs. (71) and (72) into (73) and equating coefficients of 
(i=l ,2, . . . ,6) , we obtain the generalized force components 
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( 74 ) 


^3R 


2 ^ F 2x + F 2y^2^2^ 


Q 4R " F 2z *2^ 


Q 5R 2 ^ F 3x “ F 3y^3^ £ 3^ 


3y J 


Q 6R " F 3z *3^3^ 


The radiation torque is obtained by writing simply 
3 

N r = * [r.(£ i > + u-l^.t)};^. 

Inserting Eqs. (71) and (72) into (75) we obtain the components 


(75) 


N xR " ( p lz q l " F ly q 2^1^1^ "^F F 2z q 3 + F 2y q 4^2^2^ 

" '^ F 3z q 5 + F 3y q 6^3^3^ + T ^ F 2z A 2 " F 3z A 3^ 

N yR “ F lx q 2*l^l^ + ^ T F 2z q 3 + F 2x q 4^2^2^ “ ^ T F 3z q 5 " F 3x q 6^3^3^ 



_F lx q l 4 ’l (^1 ) + 2 ^ F 2 x ‘ ^ F 2y^ q 3 <|) 2^2^ + 2 ^ F 3x + ^ F 3y^ q 5 4> 3^3^ 
+ F ly £ l " F 2x + F 2y^2 + \ ^ F 3x “ F 3y^ £ 3 


It remains to determine when these radiation effects are present, i.e., 
when solar radiation impinges on the spacecraft. This coincides with the 
period when the spacecraft is illuminated by the sun. For simplicity we shall 
be concerned only with shadowing of the spacecraft by the earth. 



Let us denote the radius of the earth by r £ and the angle between R £ 
and R c by y (see Fig. 7). Once again we assume that | R^. | » |R^| . But 

from the definitions of the cross and dot products of vectors, we can write 

I Bp x Bp I ?E ?C 

sin Y = TB£T |R c l ’ C0ST = l? E l I5 c l <77) 

Because solar radiation is in the direction of R^, we observe that the 
spacecraft is in the shadow of the earth when the projection of R^ onto R^ 
is positive and when siny < r^/ |Rq| . Hence, radiation forces will not 
be present when 

cosy > 0 and |R^ x R^| < r^ ( R £ ( (78) 

Criteria (78) can be written in the explicit form 
cosx^ (cos \p cosx^j - simj> cos i sin x^) + sin x^ cosa (cos ^ sin x^ 

+ sin ip cos i cos x^) + sin X^ sin a sin ip sin i > 0 (79a) 

and 

2 2 
{sin x^ [sin ip (sin i cos a - cos i sin a cos x^) - cos rp sin a sin x^] 

p 

+ [cos sin a sin x^ cos x^ - sin ip (sin i cos x £ + cos i sin a sin x^ sin x^)] 

+ [sin ip cos i (cos x £ cos x^ + cos a sin x £ sin x^) 

+ cos ip (cos x^ sin x^ - sin x^ cos a cos x^)] <: }' 5 < r E /R^ (79b) 
where R^ = )R e | . 


Damping Forces 


The other nonconservative forces acting on the system are the internal 
damping forces. This type of damping is generally known as structural damping, 
but the forces are often modelled as viscous forces. Of course, the forces 
being internal they produce no torques. 
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Letting c. be the damping coefficient per unit length of boom i, we 
can introduce the Rayleigh dissipation function in the form (see Ref. 1) 


F = 


1 z 

2 E 
£ i=l 



/ V 

f*1 c 

hv\ 2 h«.\ 2 

1 c i 
J o 

Vat / \at ] 


dx. 


(80) 


Using Eqs. (33), Eq. (80) becomes 


F = \ [C-j (q^ + q\) + C 2 (q^ + q\) + C 3 (q* + q|)D 


(81) 


where 

C i = f £l c i*i dx i » 1=1 .2,3 (82) 

* 0 

The damping forces can be obtained from Rayleigh's dissipation, function by 
means of the formula 

Q, D = ~ , 1=1,2 6 (83) 

3q i. 

from which it follows that 

• • • 

QlO ~ C l c <l » ^2D = C l q 2 * ^3D = C 2 q 3 

(84) 

^4D = C 2 q 4 * %D = C 3 q 5 ’ ^6D = C 3 q 6 

Quite often damping is given in terms of the damping factor rather 
than in terms of the damping coefficient c^ . The relation between the two is 

= c^/p. , i-1,2,3 (85) 

where u>\ is the lowest natural frequency of boom i and p. is the mass per 
unit length of the boom. Letting the mass of boom and balloon have the 
expression 

P^) = P i + M i 6(x. -£..), i=l,2,3 (86) 

where p.. on the right side is the constant mass per unit length of boom 
introduced in Eq. (85), is the mass of the balloon, and 6(x^ - ) is 

a spatial Dirac delta function, Eq. (34) yields 
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( 87 ) 


*i 


[p 1 - + M i <5(x. - «. i )]4> 1 - dx. = p. 4>. dx/ + (t. ) = 1 

o ■* o 


£i .2 


so that Eq. (82) yields 

2 

C i = c i — { ^- = 2^.03. [1 - M.*. 2 (£.)] , i=l ,2,3 (88) 

Equations (88) can be inserted in Eqs. (84) to express the damping forces 
in terms of the damping factors ^ . 


Mathematical Solution and Computer Program 

The numerical solution of the equations of motion involves several 
steps, namely, the determination of the lowest natural frequency and mode 
for each boom, the evaluation of certain definite integrals, the inversion 
of the matrix [A], Eq. (59), and the numerical integration of Eqs. (62). 

The computer program accomplishes all of these things with the aid of several 
subroutines contained in the IBM Scientific Subroutine Package. 

The lowest natural frequency for each boom is obtained by solving for 
the first root of the associated characteristic equation numerically. A 
simple method of interval halving produces this root. All the necessary 
definite integrals are evaluated using the trapezoidal rule in conjunction 
with Romberg's extrapolation method. Two hundred intervals are used in 
integrating each function. Matrix [A] is inverted by a standard Gauss-Jordan 
numerical procedure. An approximate solution of the first-order differential 
equations of motion, Eqs. ($2), for given initial conditions, is obtained 
by a Runge-Kutta integration procedure. Evaluation is done by means of 
fourth-order Runge-Kutta formulas with the modification by Gill. Accuracy 
is tested comparing the results of the procedure by using single and double 
increments. This method is both stable and self starting. 
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Numerical Results 


The dynamics problem of the LZEEBE spacecraft has been programmed for 
digital computation and several cases of interest have been investigated. 
These cases appear to be the most significant ones and they all differ in the 
initial conditions. These initial conditions can be divided broadly into 
two types: (1) zero (or nearly zero) initial velocities relative to an 

inertial space and zero (or nearly zero) initial displacements relative to 
an orbital system of axes, and (2) zero (or nearly zero) initial displace- 
ments and velocities relative to an orbital system of axes. These cases are 
significant because they can shed some light on the possibility that the 
spacecraft might be captured in a planar gravity-gradient stabilization 
configuration. 

The cases investigated and the results are as follows: 

Case 1 . Zero initial velocities relative to an inertial Space and zero 
initial displacements relative to the orbital axes . 

Some of the parameters have the values: = 0°, x^ = 0°, ip = 0°, 

and i = 60°. The spacecraft begins v its motion in the earth's shadow. Be- 
cause the initial tendency of the body orientation is to remain fixed in an 
inertial space, the angle tends initially to increase in magnitude at a 
rate equal to the orbital velocity n. This tends to introduce a very small 
gravity torque about the z axis. At the same time, the differential gravity 
forces on the booms cause them to deflect, with most of the elastic displace- 
ments taking place in the plane of the spacecraft. As soon as the spacecraft 
emerges from the earth's shadow, the sun's radiation pressure begins to exert 
forces and torques on the spacecraft, causing the spacecraft to rotate and 
the booms to deflect. Whereas the boom deflections remain well below one 
meter, the angles become large. In particular, the angle exceeds 2 t r. 
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thus ensuring a complete rotation of the spacecraft about its own center. 
Because of the ever-changing pattern of motion, damping appears to have no 
meaningful effect on the elastic displacements. 

In view of the large rotations of the spacecraft, it appears that the 
sun's radiation torques are sufficient to ensure that all sides of the 
spacecraft are exposed to the sun in a random-like fashion. 

For the duration of the computer run, a time period equivalent to about 
20 orbits, there was no indication that the spacecraft might settle into a 
planar gravity-gradient stabilization equilibrium. 

Case 2. Zero or small initial angular displacements and Velocities with 
respect to the orbital axes, with booms 2 and 3 deformed by differential 
gravi ty and centr i fugal forces . 

A. e 1 = 0 2 = e 3 = 0, X E = 90°, X N = 180°, \p = 0°, i = 66.5°, <* z = fi. 

This case is designed to assess the effect of the sun's radiation pressure 
on the spacecraft in the initial planar gravity-gradient equilibrium. In- 
deed, in this case the sun's radiation forces are initially normal to the 
orbit plane, and hence to the spacecraft plane, and remain constant in mag- 
nitude and direction for time intervals of the order of one day. (It should 
be pointed out that the earth's oblateness causes the orbit plane to precess 
at a rate of the order of 3° per day, so that, after a while, the radiation 
pressure ceases to be normal to the orbit plane.) The radiation forces 
cause booms 2 and 3 to oscillate about the deformed equilibrium with an ampli- 
tude of the order of 0.04 meters tip deflection. Small 9-| and 0£ oscillations 
are also present. Hence, the initial planar gravity-gradient stabilization 
equilibrium is largely maintained. It follows that, when working against 
the stabilizing effect of differential gravity and centrifugal forces, the 
radiation pressure has a very small effect. 
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B. e.| = 63 = 0, e 2 = 5°, = 0°, x N = 0°, <j> . = . o° , i . = 60°, w z = a. 

The booms undergo small oscillations about the deformed equilibrium as 
in the case 2A. In addition, the spacecraft oscillates between = ±5° 
with a period of approximately 3,000 sec. The gyroscopic effect also induces 
an oscillation about the x axis with the same amplitude but with twice the 
period. Once again the radiation pressure effect on the equilibrium state 
is minimal. There appears to be a small secular rate of reduction in 

Summary and Conclusions 

The dynamical behavior of the LZEEBE spacecraft has been investigated 
under certain simplifying assumptions. In particular, it is assumed that 
the spacecraft center of mass moves in a circular orbit around the earth, 
so that its motion in an inertial space is known. The spacecraft is subjected 
to solar radiation forces and differential -gravity forces. Thermal bending 
effects have been ignored, an assumption justified when the spacecraft 
rotates about its own mass center in a way that no one side is exposed con- 
tinuously to the sun. The formulation consists of the system Lagrange's equa- 
tions of motion for the three rotations e - ( t ) (j=l ,2,3) of the spacecraft as 
a whole and for the six elastic displacements v. (x. ,t) and w. (x. ,t) (i=l ,2,3) 

I ! II t 

of the booms, where v^ and w^ are the displacements of boom i in the plane 
of the spacecraft and normal to the plane, respectively. Note that the 
rotations 9 . are measured relative to an orbiting set of axes abc, where a 

J 

coincides with the local vertical, b is tangent to the orbit and in the direction 
of the orbital motion, and c is normal to the orbit. 1 The rotations e. define 

J 

the orientation in space of the spacecraft body axes xyz; the elastic dis- 
placements are measured relative to axes xyz. 
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FIGURE 1. LONG-TERM ZONAL EARTH ENERGY BUDGET EXPERIMENT 









The six elastic displacements are represented by one degree of freedom 

. * • ; - ♦ > . 

each, so that the spacecraft is simulated by a nine-degree-of-freedom 

* *>.■ •> . ./ 

system. For the purpose of numerical integration, the nine second-order * 

Lagrangian equations have been transformed into eighteen first-order equations 
for the state variables, namely, the spacecraft generalized displacements 
and velocities. The first-order equations have been integrated numerically 
by a Runge-Kutta procedure. Note that the differential equations are highly 
nonlinear, so that no closed-form solution is possible. 

The behavior of the spacecraft has been investigated for various cases, 
depending on the initial conditions. Computer results indicate that if the 

spacecraft is injected into orbit with zero initial displacements and velocities 

V ■*» 

relative to an inertial space, then the sun's radiation and gravitational 
torques impart to the spacecraft a rotational motion about its own center that 
can be regarded as random. The results over a 20 orbit time interval do not 
show any tendency of the spacecraft to settle into an equil ibrium corresponding •« 
to planar gravity-gradient stabilization. On the other hand, if the spacecraft 
is captured in the above equilibrium configuration, then it shows no tendency 
to escape this equilibrium state, but continues to oscillate slowly about 
that equilibrium. Hence, there is the possibility of capture in that equili- 
brium state. Note that the boom flexibility does not change the nominal con- 
figuration of the spacecraft significantly, and the thermal effects are not 
likely to change it appreciably either. Although a circular ortl It was assumed, 
it should be pointed out that orbit eccentricity has a destabilizing influence 

on the gravity-gradient equilibrium state. 

If the possibility of capture in planar gravity-gradient stabilization 

Hi 

equilibrium is to be absolutely prevented, then a slightly different spate- ^ ‘ 
craft design may deserve consideration. Indeed, a spacecraft with four equal 
booms instead of three, so that the ballons lie at the four corners of a regular 
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pyramid and the hub lies at the center of the pyr^rrtdf posfsesses suffici' 

► .4 \i * . , y ,*■: ■ 

spherical inertial symmetry to virtually eliminate gravity torques. ThV 
should permit the sun’s radiation pressure to impart to thefipacec^faft a 

rotational motion that, for all practical purposes, can be considered as 
being random. 
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